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quantum field theory on a causal set. This work can be viewed as a novel regularisation of quantum 
field theory based on a Lorentz invariant discretisation of spacetime. 
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The possibility that spacetime is fundamentally dis- 
crete presents one way to regulate the divergences in 
quantum field theory and the singularities in general rel- 
ativity. Causal set theory provides a model of discrete 
spacetime in which spacetime events are represented by 
elements of a causal set — a locally finite, partially or- 
dered set in which the partial order represents the causal 
relationships between events. 

Causal sets can be considered as physical spacetimes in 
their own right [1] (in particular, in attempting to solve 
the problem of quantum gravity [HEJIl]) or they can be 
considered as merely Lorentz invariant discretisations of 
continuum spacetimes: Lorentzian random lattices. Here 
we develop a formalism for quantum field theory on a 
causal set which is applicable to both points of view (al- 
ternative approaches include [SJH]). The main result is 
the construction of the Feynman propagator for a real 
scalar field on a causal set. This work extends previous 
results for propagators on a causal set [7]. 

The ultimate aim would be to develop physically re- 
alistic quantum field theories on a causal set. Not only 
would one expect such field theories to be better defined 
than their continuum counterparts but one would hope 
to obtain new predictions for the behaviour of matter 
which could provide evidence for spacetime discreteness. 

A causal set (or causet) is a locally finite partial order, 
i.e. a pair (C, ^) where C is a set and ^ a relation on C 
which is (i) reflexive {x ^ x); (ii) antisymmetric (x ^ y ^ 
X =^ X = y); (iii) transitive {x ^ y ^ z =4» x ^ z); 
and (iv) locally finite {\{z' € C\x < z' < y}\ < oo) for all 
x,y,z G C. Here \A\ denotes the cardinality of a set A. 
We write x ^ y if x ^ y and x ^ y. 

The set C represents the set of spacetime events and 
the partial order ^ represents the causal order between 
pairs of events. If x ^ y we say "a; is to the causal 
past of y" . The causal relation of a Lorentzian mani- 
fold (without closed causal curves) satisfies Conditions 
(i)-(iii). It is Condition (iv) that enforces spacetime 
discreteness — each causal interval contains only a finite 
number of events. 

A sprinkling is a way to generate a causal set from 
a d-dimensional Lorentzian manifold (M, g) . Points are 
placed at random within M using a Poisson process (with 
density p) so the expected number of points in a region 
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of d-volume V is pV . This generates a causal set whose 
elements are the sprinkled points and whose partial order 
relation is "read off" from the manifold's causal relation 
restricted to the sprinkled points. 

Here we shall only consider causal sets generated by 
sprinkling into d-dimensional Minkowski spacetime, M''. 
Such causal sets provide a discretisation of M"* which, 
unlike a regular lattice, is Lorentz invariant [U Sec 1.5]. 

A link between u,v £ C (written u -<*v) is a relation 
u < V such that there exists no w G C with u < w < v. 

A linear extension of the causal set (C, ^) is a total 
order (C, <) which is consistent with the partial order. 
This means that u < v =4> u < v for all u,v ^ C. 

Labelling the elements of a finite causal set (with p ele- 
ments) wi, . . . , Wp there are two pxp adjacency matrices: 
the causal matrix C and link matrix L defined by 

fj _ { I ii Vx < Vy ^ _ f 1 if <■¥ Vy 

' 1^ otherwise, ' [0 otherwise. ^ ' 



I. FREE SCALAR QUANTUM FIELD THEORY 
IN THE CONTINUUM 

We briefly review the quantum field theory for a free 
real scalar field on M'' (following the sign conventions of 
[%\). The defining equation for a Green's function G^''-' 
for the Klein-Gordon equation is 

{U + m'^)G'^'^\x)^5'^{x). (2) 

Here x = (a;°, x) = (a;°, x^) (for i = 1, . . . , d — 1), m is the 
mass of the field, 5'^ is the d-dimensional delta function 
and we choose units with = c = 1 . The d' Alembertian 
is □ := - dl 

Three important Green's functions (or propagators) 
are the retarded Gr, advanced Ga and Feynman Gp. 
In we have [71, p23]: 

Gg^(x) 0(x")0(,s2)ijo(ms), (3) 
g'^J\x) Gf{-x), (4) 
G^^\x) ^i?^')(m.), (5) 

where s = y/ (x°)2 — (x^)^ for (a;°)^ > (x^)^ and s ~ 
— jA/(a;i)2 — (a;°)2 for (x^)^ < (x^)^. Jq is a Bessel func- 

(2) 

tion of the first kind, Hq ' is a Hankel function of the 



2 



second kind and 

0(a) = 
^ we have [8, Appendix 1] 



1 if a > 
if a < 0. 



(6) 



In 



(7) 

(8) 
(9) 

where s = \/ (a;°)^ — for > if-^ and s — 

—iy/x^ — [x^Y foi' (2^'^)^ < "^1 is ^ Bessel function 
(2) 

of the first kind and is a Hankel function of the 

second kind and 5{s'^) is the delta function of . 
The PauH- Jordan function is defined as [H] 



A(a;) Gr{x)^Ga{x). 



(10) 



We drop the dimension-superscripts whenever we refer to 
arbitrary dimension. 



A. Fields in the continuum 

A free real bosonic scalar field is represented by an 
algebra of field operators (t>{x) (which act on a Fock space 
F) satisfying the following conditions: 



1. (□ + to2)0(x) = 0, 

2. ^(x) - ^{x)\ 

3. [0(a;)>(y)] = iA(y-x), 



(11) 
(12) 

(13) 



In addition there exists a Poincare invariant vacuum state 
|0) e F. With the fields so defined the Feynman prop- 
agator is given by the vacuum expectation value of the 
time-ordered product of two field operators: 



GF{y-x) = i{Q\T^{x)^{ym. 



(14) 



The time-ordering has time increasing from right to left. 
Applying (Dx + m?) to the commutator in (13) gives 



[(□, -t- m^)^{x),^{y)] = (□, + m2)zA(tj - x) = 0. (15) 



That is, even if only (12 1 and (13) hold, we have that 
+ m'^)(j){x) commutes with all the (j){y) [IT]. 



II. FREE SCALAR QUANTUM FIELD THEORY 
ON A CAUSAL SET 



Let (C, ^) be a finite causal set with p elements 
Vi, . . . ,Vp generated by sprinkling (with density p) into a 
finite causal interval in M''. 

For sprinklings into the p x p matrix 



K^^^ := aC{I-abCy 



a=l,b=- — , (16) 
2 p 



(where I is the identity matrix, C the causal matrix for 
the causal set and m is the mass of the field) is the ap- 
propriate definition of the retarded propagator [7]. 
For sprinklings into the p x p matrix 



K 



(4) 



aL(I — ahL) 



27rV6 



6 = 



(17) 



(where L is the link matrix for the causal set) is the 
appropriate definition of the retarded propagator. Both 
these matrices were defined and studied in detail in [7|. 

Our goal now is to define a, p x p matrix Kp to serve 
as the Feynman propagator on a causal set. We drop the 
dimension-superscripts to be able to refer to sprinklings 
in either or M"*. 

If Kn is the retarded propagator it follows that its 
transpose Ka ■— is the advanced propagator and 
the real matrix defined by 



A := Kr - Ka, 



(18) 



is the causal set analogue of the Pauli- Jordan function. 

The matrix iA is skew-symmetric and Hermitian. 
These two properties ensure its rank is even |10j and its 
non-zero eigenvalues appear in real positive and nega- 
tive pairs. In particular (if its rank is non-zero) there 
exist non-zero eigenvalues and normalised eigenvectors 
such that: 



iAu, 



lAu,- 



(19) 



(with \i > 0) for i = 1, . . . , s where 2s is the rank of iA. 

These eigenvectors are each defined up to a phase fac- 
tor and can be chosen such that Ui ~ v*, u\uj = v\vj — 
6ij, ulvj = for all i,j = l,...,s (where z* denotes 
complex conjugate of z and := (w*)^ is the Hermitian 
conjugate of a column vector u) . 

It's useful to define the Hermitian p x p matrix 



g := ^ A, 



(20) 



such that iA — Q 



Fields on a causal set 



We now define an algebra of field operators 4>x (acting 
on some Hilbert space i?) to represent a free real bosonic 
scalar field on the causal set. For each causal set element 
Vx (a; = 1, . . .p), we suppose there exists a field operator 
(br such that 



1. 4)^ = 



(21) 
(22) 



3. iAw = 



wx4x, = 0, (23) 

x' = l 



3 



for x,y = 1, . . . ,p. The first two conditions are natural 
generalisations of the continuum case. The last condi- 
tion ensures that any linear combination of field opera- 
tors that commutes with all the field operators must be 
zero. By (151 this is the analogue of imposing the Klein- 



Gordon equation on the field operators. 

From these field operators we can define new operators 



p 

E 

x=l 



P 

«1 '■= ^{'U'i)x(t). 
x=l 



X • 



for i — 1, . . . , s. They satisfy (aj) ' = at and 
[aijflj] = vfiAvj = —Xjulvj ^ 0, 



u, iAu.i — Xiv'ui — 0, 



(24) 



(25) 
(26) 
(27) 



and can be interpreted as (unnormalised) creation and 
annihilation operators. 

The transformation can be inverted to give 



(28) 



(here it is important we imposed (23)). 

We now define a vacuum state vector |0) € H hy 
the conditions that ai\0) = for all i — and 
(0|0) = 1. This allows us to recognise that H is the Fock 
space spanned by basis vectors {a\)"^ {a\)"^ ■ ■ ■ (aj)"=|0) 
for integers > 0, « = 1, . . . , s. 

The two-point function can be evaluated as 

s s 

{0\$.^y\0) = 5]^(«.).(^',)v(0|a^a]|0) (29) 

i=l 3 = 1 

s s 

= '^'^Mx{Vj)yXj6,j ^ Qxy. 

i=l 3 = 1 

We can now define the Feynman propagator by analogy 
with (|l4|) . With time increasing from right to left we have 



{Kp)xy '■ — i{Q\T(j)x4'y\0) '■ — * {-A-xyQyx + AyxQxy + SxyQxy) 

where A is the causal matrix for a linear extension of 
the causal set and Sxy is the Kronecker delta. In gen- 
eral there are multiple different linear extensions which 
assign an arbitrary order to pairs of unrelated elements. 
This arbitrariness does not affect Kp because the field 
operators for unrelated elements commute. 
Observing that Axy{iAxy) = {iKujxy we have 

AxyQyx — Axy{Qxy ^Axy^ — AxyQ xy ^{-^ xy ■ (31) 

(32) 



Substituting this into ( 30 1 gives an alternative form 

Kf^Kr + iQ, 
since Axy + Ayx + Sxy ^ I for all x, y = 1, . . . ,p. 



Since iA = Q — Q* the imaginary part of Q is 3(Q) = 
A/2. Combining this with (32) and looking at the real 



and imaginary parts of Kp gives 



5R(Q). 



A Kr + Ka 



(33) 
(34) 



III. COMPARISON WITH THE CONTINUUM 

The causal set propagators depend on the particular 
random causal set that is sprinkled. By calculating their 
average value for different sprinklings we can compare 
the causal set and continuum propagators. To do this, 
first fix a finite causal interval A C M'^. Pick two points 
X, y e A. Sprinkle a finite causal set into A with density 
p. Almost surely this will not contain X and Y so add 
X and Y to it to obtain a finite causal set (C, <). For 
definiteness label the causal set element X as Vi and Y 
as V2- 

We now calculate Kr and Kp for {C,<) and look at 
{Kr)i_2 and {Kf)i2, i-e. the propagator values for the 
pair {X, Y). Let ¥.{K^^^ \X, Y, W^, p) denote the expected 

value of (,K^R^)i2 (and E(A>|X, F, M'*, p) the expected 
value of {Kf)i2) averaged over all causal sets sprinkled 
into A C M'' (with X and Y added in the manner de- 
scribed and for a fixed density p). It was shown in [7] 
that 



¥.{K^^\X,Y,] 

'(4), 



p) = Gf{Y-X), 
X). 



lim E{K^,:>\X,Y,M'^,p) 



Gf{Y 



(35) 
(36) 



Using these and ( 33 1 we have 



^{E{Kf\X,Y,M^,p)) 
lim ^{E{Kf\X,Y,M'^,p)) 

p — ^OO 



^iGf{Y~X)),{37) 
^{G^p\y ~X)).{38) 



That is, the real part of the expected value of Kp is 
correct for and correct in the infinite density limit 
for [12]. We can compare the imaginary parts of Kp 
and Gp through numerical simulations. 
J By using a computer to perform sprinkling into finite 
regions of 1+1 and 3+1 Minkowski spacetime we can 
compute Kr, Ka, iA, Q and Kp for a range of sprinkling 
densities and field masses. Plotting the real and imagi- 
nary parts of Kp against the absolute value of the proper 
time for causally related and spacelike separated pairs of 
sprinkled points allows us to see if Kp agrees with Gp. 

Simulations using standard linear algebra software on 
a desktop computer give clear results for the 1+1 di- 
mensional case for causal sets with as few as 600 ele- 
ments (see Fig.[Tl. The agreement is very good provided 
<C TO <C y/p. There is disagreement between the imag- 

(2) 

inary parts of Kp and Gp as we take the field mass to 
zero but this seems to be related to the lack of a mass- 

f 2") 

less limit of ^{G'p'). There are also small deviations due 
to the calculation being sprinkled into a finite region of 
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Causally related pairs of events (p = 600, m = 10) 



Imaginary part of (Kp)^^ 
Imaginary part of Gp>(X) 




Absolute value of proper time for v^-v^ and X 



Spacelike pairs of events (p = 600, m = 10) 



Imaginary part of (K^)^ 
Imaginary part of G(?)(X) 




Absolute value of proper time for v^-v^ and X 



FIG. 1: Plots of Q{(KF)xy) for pairs of sprinkled points Vx 
and Vy in a single 600 element causal set generated by sprin- 
kling into a unit-volume causal interval in and 9(G^^ (^)) 
for X G 



Minkowski spacetime. These "edge effects" get smaller if 
the results are only plotted for pairs of points away from 



the edges of the region. 

In 3-1-1 dimensions the simulations are less clear. We 
expect good agreement with the continuum for large p 
(since the infinite density limit is needed in (36l) but 
larger sprinklings are needed to investigate this. 



IV. DISCUSSION 

The formalism provokes interesting questions about 
the eigenvectors Ui,Vi. They seem to be related to linear 
combinations of positive and negative frequency plane 
wave solutions of the Klein-Gordon equation. Plotting 
the eigenvector values as a function of position for causal 
sets sprinkled into shows relatively smooth oscilla- 
tory behaviour for the eigenvectors associated with large 
eigenvalues. They show some dependence on the shape 
of the region into which the points are placed but this 
dependence is lost when XiUiul is summed in the Q and 
Kp matrices (details to appear elsewhere). 

Extending Kp to infinite causal sets as well as further 
investigating the behaviour of Kp for sprinklings into 
M''' as well as other conformally flat or curved space- 
times would be of significant interest. Investigating the 
variance of Kp for different sprinkling densities is also im- 
portant. The propagator Kji (and therefore Kp) can be 
extended to non-sprinkled causal sets Sec 4.2]. This 
ensures both Kr and Kp are relevant to a fundamen- 
tal theory based only on causal sets with no reference to 
Lorentzian manifolds or sprinklings. 

The work presented here opens the door to interacting 
scalar field theories, scalar field phenomenology, the eval- 
uation of Feynman diagrams, S-matrices and scattering 
amplitudes on a causal set. Extending the quantum field 
theory to include interacting spinor and vector fields re- 
mains to be done — advances here could lead to physical 
predictions for matter on a causal set. 

The author would like to thank Fay Dowker for on- 
going support and encouragement. Also Rafael Sorkin, 
Johan Noldus for helpful early discussions. This work 
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Institute for Theoretical Physics. 
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